Quantum interference in nested d-wave superconductors: a real-space perspective 
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We study the local density of states around potential scatterers in d-wave superconductors, and 
show that quantum interference between impurity states is not negligible for experimentally relevant 
impurity concentrations. The two impurity model is used as a paradigm to understand these effects 
analytically and in interpreting numerical solutions of the Bogoliubov-de Gennes equations on fully 
disordered systems. We focus primarily on the globally particle-hole symmetric model which has 
been the subject of considerable controversy, and give evidence that a zero-energy delta function 
exists in the DOS. The anomalous spectral weight at zero energy is seen to arise from resonant 
impurity states belonging to a particular sublattice, exactly as in the 2-impurity version of this 
model. We discuss the implications of these findings for realistic models of the cuprates. 



I. INTRODUCTION 

Improvements in high-resolution scanning tunneling 
microscopy (STM) applied to superconductors 0, 0, 0, 
01 IE IE Ll IE E3 have raised the prospect of obtain- 
ing completely new kinds of local information about the 
cuprate materials, which may bear on the origins of the 
high-temperature superconductivity itself. Interpreta- 
tion of these experiments is understood to be a delicate 
matter, but until now has been undertaken at only the 
naivest levels for want of theoretical tools for studying 
the local properties of strongly correlated systems. As 
an example, one may consider the discovery of subgap 
impurity resonances at low temperatures in the super- 
conducting state by STMP, H Q : while comparisons 
of STM data on disordered BSCCO-2212 with the sim- 
plest calculations of a single potential scatterer in a d- 
wave superconductor 0, Qj]> E3 were understood early 
on to be only approximately successful, it was immedi- 
ately proposedjTE 0, ^| that more complicated (but 
still local) 1-impurity Hamiltonians or STM tunneling 
matrix elements could resolve the discrepancies. Only 
recently has it been pointed out that quantum interfer- 
ence of impurity states might make it difficult to observe 
true 1-impurity properties at all|l6l Il7| . In order for 
STM to fulfill its promise, it is vital to understand the 
extent to which long-range quantum interference due to 
disorder influences ostensibly local properties. 

The problem of low-energy d-wave quasiparticle excita- 
tions in the cuprates in the presence of disorder is still un- 
solved (for a review, see [Tg)- Traditionally, it has been 
assumed that the appropriate disorder potential is some 
random distribution of short-range (and possibly mag- 
netic) scatterers. More recently, there has been a grad- 
ual recognition that nanoscale spatial inhomogencities 
are frequently, and possibly always, present in HTSC. 
|E IE U El In most current theories, disorder is treated 
in the so-called self-consistent T-matrix approximation 
(SCTMA) which makes predictions for macroscopic prop- 
erties of disordered systems. The SCTMA predicts, for 



example, a constant residual Fermi level density of quasi- 
particle states p(0), which should dominate the low- 
energy transport over an energy range 7 referred to as 
the "impurity bandwidth" , in analogy to similar phenom- 
ena in semiconductors. Transport and thermodynamic 
measurements on the cuprates appear to support quali- 
tatively the predictions of this simple approach though 
there are linge ring quantititative differences which re- 
quire resolution [T3I . The SCTMA neglects "crossing di- 
agrams" corresponding to self-retracing scattering paths 
in real space, and attempts to go beyond the SCTMA 
have produced a variety of strongly model-dependent re- 
sults for the density of states (DOS), many of which do 
not support the idea of an impurity band (constant DOS 
energy range) at all. In these nonperturbative calcula- 
tions, the asymptotic limit p(0) may vanish po l pll |22| . 
saturate at a finite value [23 or diverge 0, |25l l26l de- 
pending on the symmetry of the Hamiltonian|27l |28| . We 
also note a recent semi-classical treatment of extended 
impurities suggesting a divergent density of states at the 
Fermi level |29|. 

In this paper, we perform simple, exact calculations of 
the interference of two impurities in a d-wave supercon- 
ductor, and compare to numerical calculations for many- 
impurity systems, in order to investigate the formation 
of the impurity band. Spatial fluctuations in the local 
DOS, which become quite complicated as a result of in- 
terference between impurities, contain information about 
both the SCTMA impurity band and about the quan- 
tum interference processes responsible for weak localiza- 
tion physics. For purposes of this paper, it is useful to 
make a distinction between quantum interference asso- 
ciated with weak localization and local interference pat- 
terns seen, for example, in STM experiments. We restrict 
ourselves in this initial work to a half-filled, tight-binding 
band with infinite potential scatterers. This model has 
nesting symmetries which distinguish it from the cuprate 
superconductors, but is nevertheless interesting from two 
points of view, that transparent analytical results for 
some properties can be obtained, and that the charac- 
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ter of the divergence of the density of states near half- 
filling is controversial The two- impurity problem is 
the simplest problem which includes the interference pro- 
cesses which lead to the formation of the impurity band, 
as well as processes which lead to weak localization. 

Early work on the two-impurity problem in a d-wave 
superconductor was numerical in nature and focussed on 
the local density of states (LDOS), exhibiting unusual 
local interference patterns which depended on the orien- 
tation of the vector R separating the two impurities. [30l | 
More recently, the relation to impurity band formation 
was discuss ed l3ll a nd predictions were made for STM 
experiments [la. fl7| . assuming that "sufficiently isolated" 
two im pur ity configurations could be identified. In ref- 
erence pjj, the bound state wavefunctions of the two- 
impurity system were identified and classified. By anal- 
ogy with the molecule problem in quantum mechanics, 
one expects that the single-impurity resonance energies 
split as the impurities are brought together, and that the 
wavefunctions are formed from symmetric and antisym- 
metric combinations of the isolated impurity wavefunc- 
tions. In fact, because of the particle- hole and fourfold 
rotational symmetries of the superconducting state, the 
situation is more complicated, with the effective overlap 
depending on R . Indeed, it has been shown that for 
many pair configurations, the density of states does not 
consist of four well-defined resonances pH Il7| . 

The interference between impurities persists up to 
large impurity separations. In Ref. |17| it was noted that 
two impurities with R || (110) could cause splittings com- 
parable to the original resonance energy for R of many 
tens of lattice spacings. The spatial LDOS maps are 
therefore very different from superposed single-impurity 
maps, and one may ask the question whether this distinc- 
tion persists in the case of many impurities. That is, is it 
to be expected at experimental impurity concentrations 
that a resonance found by STM really corresponds to an 
isolated impurity whose LDOS is predictable within a 
simple 1-impurity model 17]? Alternatively, arc interfer- 
ence effects omnipresent, destroying expected 1-impurity 
resonances and leading to new, long range LDOS pat- 
terns which require a many-impurity interpretation? If 
the latter scenario is realized, how can it be that STM 
experiments seem to see such similar spectra on or near 
impurities embedded in very different local disorder en- 
vironments? We resolve these questions below by argu- 
ing that in the generic case the individual many-impurity 
eigenstates are highly distorted from mere superpositions 
of 1-impurity LDOS patterns, but that STM measure- 
ments tend to average over many such eigenstates, can- 
celling some of the long-range effects of interference. Ex- 
ceptions are very low energy states of the nested d-wave 
superconductor, which experience symmetry-driven level 
repulsion effects which prevent such cancellations. These 
considerations lead to a picture where, with the excep- 
tion of the zero energy states, the local impurity reso- 
nances appear homogeneously broadened to any probe 
which averages over a macroscopic energy window. This 
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FIG. 1: Schematic figure of the many-impurity DOS (a) in 
the unitary limit of the half-filled band and (b) in the unitary 
limit of a generic band. The plateau in the impurity band is 
characterized by a nearly constant density of states po. The 
zero-energy suppression in (b) is discussed in Refs. |2l|, H3 



result has important consequences for the interpretation 
of STM spectra. It means that, while the position of a 
peak and the crude local LDOS pattern at an energy near 
the peak may indeed qualitatively reflect one-impurity 
properties, e.g. the strength of the 1-impurity potential 
Vo, the widths of spectral features measured at any site 
will reflect the impurity bandwidth 7 characteristic of the 
disordered system as a whole. 

The second goal of this paper is to investigate the di- 
vergence in the total density of states in the completely 
nested model from a local point of view, applying what 
we can learn about the 2-impurity system. In contrast to 
the prediction of a residual p(0) by the SCTMA, Pepin 
and Lee (PL) [2^] found that, for an infinite scattering 
potential, the disorder-averaged density of states should 
diverge at the Fermi level (taken to be the zero of en- 
ergy here) as p(uj) « rii /uj log 2 w where ni is the density 
of impurities. The schematic picture of the total DOS 
in this case is exhibited in Fig. la), to be contrasted 
with the more generic case expected in the absence of 
the nesting symmetry (lb). One surprising aspect of the 
PL result for the TV-impurity T-matrix is that it is es- 
sentially rii times the single- impurity result. Upon closer 
inspection, however, their result is not directly tied to 
the one-impurity resonance at the Fermi-level, but is the 
result of interference between distant impurities. Nu- 
merical calculations showed that the divergence arises 
because of a global particle-hole symmetry [23 which is 
particular to the tight-binding model at half-filling. It 
was later shown that this nesting leads to a novel diffu- 
sion mode j2^| producing a positive logarithmic correction 
to the DOS. This general structure of the divergence has 
also been found by Mudry and Chamon|25| and numeri- 
cal calculations 32J seem to confirm it, although in both 
cases the strength of the divergence could not be veri- 
fied. The situation is not settled, however, and other re- 
cent field-theoretical approaches 26J find a different form 
for the divergence which is reminiscent of the half-filled 
normal metal. The investigations of the 2- and many- 
impurity problem presented here paint yet a different 
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picture. While (for reasons discussed in the text) it is 
difficult to rule out the existence of a continuous diver- 
gent contribution, we argue that the strong divergence in 
p(E — > 0) seen in previous numerical work is actually in- 
dicative of a delta-function divergence at the Fermi level. 
In the final stages of writing this work, we became aware 
of a recent conserving weak-localization calculation!^ 
which comes to the same conclusion. 

The paper is organized as follows: in section lll Al we de- 
rive expressions for the Green's functions G(r, u>) needed 
to evaluate the two-impurity T-matrix. Asymptotic ex- 
pressions for large r have been found previouslyj^ Ell 
l35| and our results, valid for small u>, are complementary. 
In section III Bl we specialize to the fully nested strong 
scattering model, evaluate the density of states for dif- 
ferent impurity configurations, and show that there are 
three different classes of impurity-pair orientation. For 
two of these classes, the DOS diverges as u> — ► 0, while 
the DOS vanishes for the third. In all three cases, inter- 
ference between impurities is substantial as in Ref. |24j . 
but ultimately the observed divergences arise from the 
local rather than nonlocal correlations. In section ITTT1 we 
establish a connection between the zero-energy LDOS of 
the fully nested disordered system and the zero-energy 
DOS of the one- and two-impurity problems. For the 
fully nested model, we find that, in a given configura- 
tion, only impurities on a given sublattice contribute to 
the resonant weight at zero energy. The impurities in this 
class form a network with spatial separations equivalent 
to the resonant configurations in the 2-impurity case, and 
numerical scaling of the total spatially integrated DOS is 
shown to be consistent with p{uS) ~ 

In section ITVl we summarize our conclusions and dis- 
cuss the less symmetric situation found in the cuprates. 
We argue that, because the STM averages over many 
multi-impurity eigenstates, the LDOS indeed appears to 
represent a set of nearly isolated impurity states with 
spectral features which are similar from impurity to im- 
purity. On the other hand, we expect the width of these 
local states in energy to be typically the impurity band- 
width arising from the full disordered system. 



II. TWO IMPURITIES IN A HALF-FILLED 
BAND 

A. Green's functions 

The BCS Hamiltonian for a pure <i-wave singlet su- 
perconductor in a tight-binding band can be written as: 

Ho = 53$£[(ek-A»)T3+A fc fi]$ k , ( la ) 



£k = — 2t(cos k x + COS fcy), 
Ak = Ao(cosfca; — cosfc y ), 



(lb) 

(1c) 



where $k = ( c kj.c^ k |) is a Nambu spinor, and ft are the 
Pauli matrices. Energies are measured relative to the 



center of the band, so a chemical potential of p = cor- 
responds to half-filling. The associated Green's function 
is, in real space, a function of the relative coordinate 
r = (m,n), where r is measured in units of the lattice 
constant and m and n are integers: 



y e ik-iv-,n 



G°(k,w), 



cos(k x m) cos(k y n) 



utq + (efc - p)r 3 + A k fi 



(2) 



where £k = \/ e k + ^k denote quasiparticle energies, 
denotes a matrix in Nambu space and the superscript 
denotes the bare Green's function. Frequently, it is con- 
venient to make the decomposition in terms of Nambu 
spinors 

0=0 

An ensemble of N short-range scattering potentials at a 
set of sites R^ introduce a perturbation 



Hi, 



V 



N 



i=l 



R 



f 3 $ 



R, 



where Vo is the strength of the impurity potential. For- 
mally, there is an exact solution for the disordered Green 
function in terms of the 2N x 27V many-impurity T- 
matrix: 



G(r,r',w) = G°(r-r',w) 

]G°(r - Ri^fyM&iRj - T,u) 



with i,j the position indices of the impurity sites, and 

f = [1 <g> f - fgVoGV)]- 1 ! ® hVo, 

where the boldface indicates a matrix in spatial indices 
in the subspace of impurity sites (ie. G°(w) = G°(Ri — 
Rj,cj)) and the inverse is a matrix inverse. In the limit 
of a single impurity, the T-matrix simplifies to T(u>) — 
[V^fa-G ^,^)}- 1 , withG°(0,w) = G°(r = 0,w). This 
limit has been studied extensively. 

In this work, we are particularly interested in the two- 
impurity T-matrix with one impurity at the origin (for 
simplicity) and the other a displacement R = (m, n) from 
the origin. The 2-impurity T-matrix is a 4 x 4 matrix 
which satisfies 



V^3-G (0,^) -G°(R,«) 
-G°{TL,oj) V^3-G (0,w) 



(3) 



Expressions for the local Green's function G°(0, ui) have 
been derived in many places, but the nonlocal Green's 
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function G°(R, uj) is less well understood, although sev- 
eral asymptotic expressions have been found [24l 134 135 | . 
In the Appendix, we derive expressions which are valid 
for the half-filled band, and which become exact in the 
limit uj — > 0. 

We find that the local Green's function for general 
complex uj is 



G°(0,w) 



aoj , A 2 
— — In 



:To, 



(4) 



where a = N/(2ttvfva), N = 4 is the number of nodes, 
vf is the Fermi velocity and va is the anomalous quasi- 
particle velocity |VkAk|, and the cutoff A is of order A . 
The expansion in uj for r = (to, n) depends on whether 
n and to are odd or even. For the (even, even) case, we 
have 

G°(r,w) ^ (-1)^ [G (0,uj)+ujC (r)]f , (5) 

where Go(r) is a real function of r. We find similar 
leading-order expressions for (m,n) — (odd, odd), 

G Q (m,n,uj)^ujC Q (r)T , (6) 

while for (m,n) = (odd, even) or (even, odd), 

G°(m,n,uj) -» Ci(r)fi + G 3 (r)f 3 , (7) 

where Gi(r), and G 3 (r) are real constants. This distinc- 
tion between even and odd sites accounts for the oscilla- 
tory nature of the wavefunctions for the special case that 
the Fermi wavevector is commensurate with the lattice. 



B. Density of states for two impurities 

In this section, we derive expressions for the density of 
states for two impurities in a half-filled band. The discus- 
sion focusses on the unitary limit Vq — ► ±oo. The half- 
filled tight-binding band possesses a particular global 
nesting symmetry [H| T2G°(k + Q, uj)t2 — — G°(k,u;), 
with Q = (tt, 7r) which is satisfied at uj — 0. For simplic- 
ity, we call this the Ti symmetry. Potential scattering 
violates this symmetry, but in the case of infinite po- 
tential, impurity sites are effectively removed from the 
lattice, and the symmetry is recovered for any disorder 
configuration at uj = 0. In real space (see eg. 26])the r 2 
symmetry may be expressed as 



T 2 G(r,r',uj)T 2 



_ _ iQ-(r-r') 



G(r,r',uj). 



(8) 



It will be useful to decompose the square tight-binding 
lattice into the usual two interleaved sublattices (denoted 
A and B). The phase factor on the right hand side of Eq. 
is +1 if r and r' belong to the same sublattice, and 
— 1 otherwise. 

The simplest quantity of interest is the quasiparticle 
density of states 

p(uj) = - En) + 8(uj + E n )] 

n 

= PoM + Sp(uj), 
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FIG. 2: Change in the quasiparticle density of states arising 
from impurities separated by R = (2, 2) as a function of en- 
ergy w/t for scattering potential Vo = 10 6 t. Inset: Scaling of 
the resonance peak energies as a function of Vo. The DOS for 
R = (2, 0) is almost identical. Energies are measured in units 
oit, and A = O.lt 



where E n are the positive energy eigenvalues of the su- 
perconducting Hamiltonian, po(uj) is the DOS of the 
disorder-free system and Sp(ui) is the change induced by 
the impurities. The DOS is related to the 2-impurity 
T-matrix defined in Eq. @ by the phase shift 77(0;) |3iJ: 



Sp(uj) 



1 drj 

7T 8uJ ' 



(9) 



Im dctT 



(10) 



where r\ is given by, 

ri(uj) = tan" 

Re detT 

and the determinant is over spatial and spin indices. 

We start with a discussion of two impurities belonging 
to one of the sublattices. The two impurities are at R4 
and R2 with R = Ri — R2 = (to, n) =(even,even) or 
(odd, odd). The two- impurity T-matrix defined in Eq. 
is particularly simple in this case: 



T = I 
D 



-Gg(0,w)fb G° (R,uj)t 
Gg(R,cj)f -G8(0,w)fb 



where D = G[j(0,w) 2 - Gg(R, uj) 2 . Noting that 
detf = ±, 

we keep the leading order terms in Go(R, uj) as uj — > 0, 
given explicitly in Eq. Q and JHJ and find that det T 
diverges as 

f [2wG (R)G[|(0,w)]- 2 R = (even, even) 
\G° (0,uj)- 4 R= (odd, odd) 

and (analytically continuing uj to the real axis) 



detT 



6p(uj) 



1/[uj log 2 (A/o;)] R = (even, even) 
2 / [uj log 2 (A/ uj)} R=(odd,odd) 



(11) 
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that det T = D' 2 and that 

Sp(u 0) oc 4- hi - | -> (12) 
aw \ w / 

A similar result holds for R = (odd, even). Physically, the 
fact that Sp vanishes at the Fermi level indicates that 
bound state energies must always arise at nonzero ener- 
gies. Numerical calculations of the DOS shown in Fig. [3| 
demonstrate that there is no remnant of the single impu- 
rity uj — > divergence for this orientation, and that the 
resonance energies scale very little with Vq. In this case, 
it is the dominance of the nonlocal terms which shifts the 
resonance to finite energy. 



FIG. 3: Change in the quasiparticle density of states for impu- 
rities separated by R = (2, 1) as a function of energy ui/t for 
scattering potential Vo = W 6 t and Ao = O.lt. Inset: Scaling 
of the resonance peak position as a function of Vq. 



Because of the similarity of the approaches, we are able 
to compare our findings with those of PL [24| in some de- 
tail. Although the form of Eq. (|11|) is suggestive of the 
asymptotic result of PL for the disorder-averaged limit, 
its origin is quite different. This difference is easiest to 
see for the (odd, odd) impurity configuration: here the lo- 
cal Green's function G (0, u>) is dominant over the non- 
local term G°(R,u>) and the physics of the low energy 
resonance is essentially that of two non-interacting im- 
purities. The total weight of the resonance is therefore 
twice that of a single impurity. For the (even, even) case 
the situation is a little more complicated, since the local 
and nonlocal terms are nearly equal in magnitude; inter- 
ference effects reduce the spectral weight of the combined 
resonance to half that of two isolated resonances. In both 
cases the situation is quite different from Ref. |2J] where 
the logarithmic divergence arises from averaging over all 
possible impurity separations using the approximate form 
G°(R, u>) ~ 1/R out to a cutoff - t/R. The PL result is 
inherently nonlocal. 

Numerical calculations for two impurities with separa- 
tion R = (2,2) are shown in Fig.0 For V a = lOOt, four 
clearly defined peaks are seen, corresponding to the level 
splitting of the single impurity resonances of the isolated 
impurities As shown in the inset, the peak positions 
scale strongly with Vo, and a single peak appears only 
when V Q ~ 10 5 t. 

We continue now with the case where the impuri- 
ties belong to different sublattices and are separated by 
R = (even, odd). The two-impurity T-matrix defined in 
Eq. is: 



III. DISORDERED SYSTEM WITH GLOBAL 
PARTICLE-HOLE SYMMETRY 

In this section, we discuss the correspondence between 
the two impurity problem and the disordered <i-wave su- 
perconductor. There are two separate issues to be dealt 
with. The first has to do with the nature of the diver- 
gence at u> = which occurs in the tight-binding model, 
while the second has to do with the more general question 
of how the impurity band evolves with impurity concen- 
tration. For these calculations, we numerically diagonal- 
ize the mean-field Hamiltonian for a random distribution 
of impurities, under the assumption of a homogeneous or- 
der parameter for a finite sized LxL system with periodic 
boundaries. For a detailed description of the method, we 
refer the reader to eg. |2^. We retain the eigenenergies 
E n and the eigenvectors 



(")fV> - 



iJrW(r) 



u(")(r) 
!)W(r) 



E„ 



The total density of states is just p(w) 
and the single-spin LDOS is 

p(r,wH ]>>("> (r)| 2 ^-i;„). 

n 

Since there is no moment formation, a =] and a =J, are 
equivalent. 

The Green's function G°(k, uj = 0) for the Hamil- 
tonian (|la|l (with p = 0) has the special symmetry 
f 2 G°(k + Q,0)f 2 = -G°(k,0) where Q = (tt.tt) is 
the antiferromagnetic wavevector. The t 2 symmetry is 
required[22l for the divergence in p(u> — > 0). This 
symmetry is only strictly satisfied when L is even [37], and 
an even-odd oscillation at the Fermi level as a function of 
L is clearly evident in our numerical work. Throughout 
this paper, we restrict ourselves to even L. 



T = — 

D' 



-Gg(0,cj)f G 1 (R)f 1 + G 3 (R)f 3 
G 1 (R)f 1 + G 3 (R)f 3 -Gg(0,cj)f 



with D' = G(0, to) 2 - Gi(R) 2 - G 3 (R) 2 . It follows easily 



A. Divergence at u = 

The DOS for a large concentration m = 0.1 of 
strong scattering impurities in a d-wave superconductor 
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FIG. 4: Total density of states for n, = 0.1. (a) DOS 
for Vo = WOt (solid) and L = 60. Eq. JHJ is plot- 
ted for comparison (dashed line), (b) Scaling of DOS with 
Vo- ppl(cj) is again plotted for comparison, (c) Scaling 
of DOS with L. (d) Scaling of DOS with Vo for Vo/t = 
100, 500, 1000, 5000, 10 4 , 10 5 , 10 6 and L = 60. A back- 
ground po = 0.25t -1 has been subtracted. The figure shows 
that the density of states is a peaked function whose width 
scales as 1/Vb and whose height scales as Vo, suggesting that 
limv-Q^oo p(uj) ~ S(ui). All energies are in units of t. 



is shown in Fig. ^ The figure is restricted to low ener- 
gies, and shows only the zero-energy peak at the Fermi 
level, and a small portion of the impurity band. For com- 
parsion, the d-wave gap has an energy Ao = 0.2t and the 
gap edge in the tunneling density of states is OAt. For 
clarity, we often make a distinction between states in the 
peak and states in the impurity band, by which we mean 
states belonging to the DOS plateau which is character- 
ized by a constant density of states pq. In Fig.^Ja), for 
example, p w 0.25i -1 . 

In Fig. Ufa), the total DOS is shown for an impurity 
potential Vq = lOOi corresponding to a strong scattering 
potential. The results are in quantitative agreement with 
earlier numerical work|27l|33- The PL result 



Mpn'(A/w) + (tt/2)2 



(13) 



is also shown. Here, we take A = 1, first because this 
was the cutoff used in previous numerical work 32J and, 
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FIG. 5: Scaling of the peak area and inverse participation 



ratio with system size L and Vo 



10°t. 



is averaged 



over states with energies 10" 5 t < E n < 0.031 (a) Scaling for 
rii — 0.1. (b) Scaling for m = 0.2. Solid lines are linear fits 
to the data. For these curves, Ao = 0.5t. 



second, because this gives a good fit to the numerics at 
Vo = lOOt. It should be clear from Figs. Ufa) and (b) 
however, that although the fit is striking at Vq = lOOi, 
it is less so for other values of Vq- In our numerics, we 
find a smooth evolution of the low energy peak as a func- 
tion of Vo and there is no value of Vq beyond which the 
asymptotic behaviour saturates. In general, ppl{u) does 
not appear to fit the data well, except for certain spe- 
cial parameter sets. The shape of the peak at lo = is 
modified by finite-size effects. There is a crossover in be- 
haviour which occurs when the mean level spacing in the 
impurity band 5l = l/(poL 2 ) is comparable to the peak 
width. Scaling of the DOS is shown for Vq = 250t in Fig. 
life). The peak height scales with L for L < 40 and satu- 
rates at larger system sizes. The implication is that some 
care must be taken in approaching the Vq — * oo limit. 

The unitary limit of the infinite system may be ap- 
proached in two ways. First, one may consider taking 
limy^oo hin^^oo so that the level spacing in the impu- 
rity band is much less than the peak width. Second, one 
may consider taking the limit L — > oo with Vq = oo. In 
the first approach, the Ti symmetry is only strictly sat- 
isfied when L = oo, while in the second approach, the T2 
symmetry is rigorously satisfied for any even value of L. 
For this reason, we view the second approach as prefer- 
able. 

The limit Vq — > oo for fixed L is illustrated in Fig. 
0fd). The data are scaled by the impurity potential, 
and the general trend is that as Vq is increased, a sharp 
peak develops at to = 0. Furthermore, the peak scales as 
p(uj) w VqF^lvVq), implying that 



lim P (uj) 

Vo— >cx> 



(14) 



Not surprisingly, the weight contained in the delta- 
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FIG. 6: Local density of states for 0.4% concentration of im- 
purities and \E„\ < 10 _5 t (4 eigenvalues). Impurity loca- 
tions on sublattice A are indicated with open circles, those on 
sublattice B with filled circles, and the impurity potential is 
Vo — 10 6 i. Inset: The inset shows a detail of the LDOS for a 
single impurity. 



FIG. 7: Local density of states for 2% concentration of im- 
purities and \E n \ < 10~ J t (20 eigenvalues). Impurity loca- 
tions on sublattice A are indicated with open circles, those on 
sublattice B with filled circles, and the impurity potential is 
Vo = 10 6 t. 



peak in the Vo — > oo limit scales with L, as shown in 
Fig. El For rii = 0.1, this scaling is consistent with what 
we found in Fig. E{c). When = 0.2, on the other 
hand, the peak area saturates when L > 40, which is 
not expected since the peak width is still many orders 
of magnitude smaller than the typical level spacing 5l 
in the impurity band. To learn more about the origin of 
this saturation we plot in the same figure the scaling of 
the inverse participation ratio, defined by 

^ EJ^Hr^W")^) 4 ] 
a(ui) = > k6[oj — E n ). 

» (Ei[« (n) ('0 2 +« (n) ('0 2 ]) 

a(uj) scales as L~ d for wavefunctions which are extended 
in (i-dimensions, and does not scale with L for localised 
states. The localization length is typically extracted from 
the crossover which occurs when L w £z, where £l is 
the localization length. As we shall see below, states in 
the delta-peak behave differently from those in the im- 
purity band, and we find that the peak area is correlated 
with the localization properties of the impurity band. In 
Fig. El the inverse participation ratio is averaged over 
states in a narrow energy window adjacent to (but not 
including) the delta peak. It is evident from the figure 
that for m — 0.2, a crossover to the localized regime oc- 
curs, and we can extract a localization length £t sa 40. 
Remarkably, we find that the area of the <5-peak appears 
to saturate when L > This situation is analogous 



to one reported earlier in c?-wave superconductors pos- 
sessing no special symmetries. There, it was shown that 
quantum interference (arising from "maximally crossed" 
diagrams) leads to a suppression of the DOS at the Fermi 
level pi) over an energy scale Sg L = l/(po£|) (This situa- 
tion is illustrated in Fig^b)). In finite size systems, the 
energy scale for the DOS suppression is actually 5l and 
the scaling of the suppression saturates when L > £l [27| • 

We note that the origins of the delta-peak divergence 
are fundamentally different from those discussed in PL, 
where the divergence arises from the cumulative effects of 
interference between a large number of distant impurities. 
Here the result appears to be a mesoscopic effect which 
survives because localization makes the effective system 
size finite. However, although the delta-peak result is dif- 
ferent from earlier predictions for a continuous divergence 
at the Fermi level, it does not preclude the existence of 
an additional divergent term which is unobservable be- 
cause of finite system size effects. Indeed, if we consider 
the effect of finite system size on the PL result we find 
the interference between distant impurities is cut off by 
L and we should make the substitution lu — * max(u;, t/L) 
in ppl, implying a cutoff energy oj c ~ t/L below which 
the DOS saturates. By this estimate, the contribution to 
the plot in Fig.^Jis cut off below lu c ps 0.017£, suggesting 
that the PL peak should be unobservable. 

It is particularly instructive to consider the structure 
of the delta-peak divergence in real space. Figure[S]shows 
the combined local density of states from the eigenstates 
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FIG. 8: Local density of states for 0.5% impurities and Vb = 
I0 6 t derived from a single eigenvalue with E„ — 0.0318£. 



with energy \E n \ < lCP 5 i which make up the delta-peak 
(these states are well-separated from all other eigenval- 
ues). For a single impurity (shown in the inset) the 
zero-energy resonance has a fourfold spatial structure 
with bright lobes on sites adjacent to the impurity along 
the antinodal (100) and (010) crystal directions, and ex- 
tended tails in the nodal (110) and (110) directions, in 
agreement with many earlier calculations [llj. For 0.4% 
disorder (10 impurities), the situation is quite different; 
even at this relatively low concentration, there is signif- 
icant interference between impurities. We see four pro- 
nounced zero-energy resonances, but the remaining six 
impurities are — at best — only weakly visible. For each 
of the visible resonances, the LDOS has the superficial 
structure of the isolated impurity LDOS, with maxima 
appearing in the antinodal direction and tails extending 
away from the impurities in the nodal directions. How- 
ever, there is no obvious correlation between the degree of 
isolation and the appearance of a zero-energy resonance. 
Indeed, of the four strong resonances, only two are more 
than 10 lattice sites from the nearest impurity. For 2% 
disorder (50 impurities) , shown in Fig. [7| the situation is 
similar. Only a small fraction of impurities contribute to 
the zero-energy LDOS and, again, the visible resonances 
do not necessarily belong to the most isolated impuri- 
ties. At this higher impurity concentration, however, a 
definite pattern in the LDOS is observable. Long tails 
along the (110) and (110) directions give the appearance 
of a network of impurities. 

Remarkably, we find that all impurities within the visi- 
ble network in Fig.0belong to one sublattice, arbitrarily 
denoted A, while the remaining impurities belong to the 



B sublattice. Similarly, in Fig. |B| all visible impurities 
belong to the B sublattice. While this is reminiscent of 
the two-impurity problem discussed in the previous sec- 
tion, it is also quite surprising. For the two-impurity 
problem, it was shown that the zero-energy resonance is 
preserved when both impurities inhabit the same sublat- 
tice and is destroyed otherwise. The natural extrapola- 
tion is that, for a random distribution of many impurities, 
every impurity is expected to have have some reasonably 
close neighbor belonging to the other sublattice which 
contributes to the destruction of the zero-energy peak. 
Clearly, this does not happen. Instead, the impurities be- 
longing to the A sublattice for this sample are dominant 
at lu — for reasons we do not completely understand 
at present. An apparent consequence of this dominance 
is that the resonances of impurities belonging to the B 
sublattice are shifted to higher energies. We speculate, 
but cannot prove, that the system in the thermodynamic 
limit will have "domains" of typical size £l in which ei- 
ther A or B impurities are resonant. 

The observed networks are also reminiscent of an ear- 
lier proposal|35j that impurities form networks from sin- 
gle impurity resonances which lead to a derealization 
transition as to — » 0. Numerical scaling calculations |32j 
for a finite impurity potential (Vb = 100£) did not find 
such a transition, however, nor does the present work 
(see below). In any case, we emphasize that the sharply 
defined networks exhibited above are a feature of Hamil- 
tonians with ti symmetry only, and not a general feature 
of d-wave superconductors as suggested in |35| . 



B. Impurity band away from u — 

We now turn our attention to the states in the impurity 
band away from lu = 0. The previous analysis raises some 
interesting questions about the formation of this "band" . 
Is the u> > DOS plateau formed, as Figures|Hland[7|pcT- 
haps suggest, by summing over many impurities, some of 
which are resonant at a given energy and others not? 
This would imply that, as energy was scanned in STM 
experiments, different impurities would "light up" - be- 
come resonant-and turn off at different energies within 
the impurity band, a scenario we will refer to as "inho- 
mogeneous broadening" of the impurity resonances. Ex- 
perimental data indicate instead that all impurities, 
regardless of local environment, appear to be resonant all 
through the impurity band, i.e. that each local spectral 
function is qualitatively similar in position and width, 
i.e. "homogeneous broadening" . In addition, there is 
some evidence from explicit Zn substitution |j| that the 
number of impurity resonances corresponds closely to the 
number of Zn atoms introduced into the crystal, i.e. there 
are no atoms which do not light up. It is for this reason 
that interpretations have typically been given in terms 
of one-impurity models. However, in the same experi- 
ments the width of spectral features is roughly an order 
of magnitude larger than those predicted by the simplest 
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FIG. 9: Local density of states for 0.5% impurities and Vb = 
10 6 t averaged over 5 eigenvalues in the energy interval \E n — 
0.03*1 < 0.02t 
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FIG. 10: Inverse participation ratio a(u>) for n, = 0.1, Vb = 
10 6 on a 30 x 30 lattice with 50 configurations. States in the 
impurity band for \E n \ > t/Vo have approximately uniform 
spatial extent. States with [i5 n |Vo < t exhibit strong fluc- 
tuations in spatial extent. Note that states with a(E n ) = 1 
are confined to a single site, while states with a(E n ) ~ iV^ 1 
have a uniform spatial distribution. Inset: A histogram of the 
distribution of a(ui) for states in the zero-energy 5-function. 
Note the logarithmic horizontal axis. 



one-impurity models. 

These apparent paradoxes can be resolved by recog- 
nizing that the energy range probed by STM, although 
very small (O(O.lmeV)) in laboratory terms, is still large 
enough to sample an essentially infinite number of eigen- 
states of the macroscopic system. In Fig. [S] we show 
the LDOS derived from a single eigenstate at an energy 
which is in the impurity band, but away from the zero- 
energy delta-peak. Two features of this figure stand out. 
First, as was the case at w = 0, only a fraction of the im- 
purities contribute to any given eigenstate. Second, the 
extended tails which were are important in the formation 
of the delta-peak are blurred by the incommensurability 
between the lattice and the wavevectors contained in the 
eigenstate. As we move further away from 10 = 0, this 
incommensurability becomes more pronounced and the 
tails become increasingly blurred. 

The inequivalency between impurities in Fig. |H1 is sur- 
prising not only because STM provides little evidence for 
such a picture, but also because the arguments about the 
formation of networks fail when u) ^ (indeed, there is 
no visible network in the figure). When one now aver- 
ages the LDOS over a small energy window, as in Fig. 

the system starts to look much more homogeneous, in 
the sense that all impurities contribute visible resonances 
with the classic fourfold symmetry. The window width 
is small compared to the impurity band, and we have 
checked that the pattern averaged in this way remains 
roughly the same up to energies of order the impurity 
band itself, 7 ~ 0.25t for the parameter set of the figure. 
Thus, it appears as if there is an important distinction 



between individual cigcnstatcs which determine, for ex- 
ample, localization properties, and averages over finite 
energy windows which determine the tunneling spectrum. 

Finally, in order to solidify the connection between lo- 
cal and bulk properties of the disordered system, the en- 
ergy dependence of the inverse participation ratio is plot- 
ted in Fig. llUI For each impurity configuration, a(E n ) is 
calculated for all the cigcnstatcs in the spectrum, and the 
aggregate is shown for 50 impurity configurations in the 
figure. There is a clear distinction between states inside 
and outside of the <5-peak. States outside the (S-peak are 
clearly extended (the localization length is much larger 
than the system size) and the distribution of a(ui) is rela- 
tively narrow at a given energy. On the other hand, there 
is a broad distribution of a(u>) in the <5-peak, indicating 
a mix of localized and extended states. The figure inset 
shows a histogram of the distribution that demonstrates 
that most of the spectral weight in the S-pe&k comes from 
the extended tails of the resonances (Fig. 0) and not from 
the highly- visible localized resonances. 



IV. CONCLUSIONS 

In this work we have studied the unitary limit of a 
disordered, half-filled d-wave superconductor with tight 
binding band. This model has a particular symmetry 
which is known to lead to a divergence in the density of 
states at the Fermi level, although the particular form 
of the divergence is controversial. We began with a dis- 
cussion of the two-impurity problem, which yields an an- 



10 



alytical solution in the uj — > limit. We found that, 
owing to the commensurability of the nodal wavevectors 
and the tight-binding lattice, there is an even-odd oscilla- 
tion in the two-impurity density of states in the unitary 
(infinite scattering potential) limit. Impurity pairs on 
the same sublattice have a zero-energy divergence in the 
DOS similar to the single-impurity divergence. The ori- 
gin of this divergence is quite different from that reported 
earlier . which arises from the cumulative interference 
of a large number of distant impurities. 

We also noted that for impurities located on differ- 
ent sublattices, the zero-energy single-impurity resonance 
is shifted to higher energies as a result of interference. 
Based on this result alone, it is natural to assume that, 
in the many impurity limit, any remnant of the single 
impurity peak will be obliterated since each impurity is 
expected to have at least one reasonably near neighbor 
which lies on the other sublattice. Surprisingly, we found 
that this is not the case. Exact numerical studies of finite- 
size systems show that unitary impurities actually form 
two interleaved networks on the A and B sublattices, 
one of which contains spectral weight at u> = 0, while 
the other does not. Away from u) = 0, quasiparticle 
eigenstates are no longer commensurate with the lattice, 
networks connecting resonant states along the nodal di- 
rections are smeared, and individual eigenstates consist 
of distorted resonances, which are inhomogeneously dis- 
tributed. When the LDOS is averaged over a small win- 
dow in energy, however, as in an STM experiment, the 
fourfold nature of the 1-impurity resonances is qualita- 
tively recovered, and resonances on individual impurity 
sites appear remarkably similar, provided the impurities 
are not in immediate proximity. Although the resonance 
peak positions may be qualitatively related to the res- 
onant energies of the underlying 1-impurity model, the 
widths are very different, of order the impurity band- 
width, given in the unitarity limit by 7 ~ ni\/ AoEp. 

The uj > states of the tight-binding band are generic 
in the sense that they do not posses the T2 symmetry, 
or any other symmetry which is not present in high 
T c superconductors. In this sense, our results should 
be qualitatively applicable to the experiments on real 
cuprate materials. They suggest that the ability of one- 
impurity models of any kind to explain the details of 
local STM spectra in samples with per cent level disor- 
der are severely limited. To substantiate this picture, 
it will be useful to compare local spectra on sites (e.g. 
impurity or nearest neighbor sites) around different im- 
purities using realistic bands. Numerical calculations to 
realize the large systems necessary to obtain the resolu- 
tion required to reach definite answers to these questions 
are in progress. 
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APPENDIX 



The purpose of this appendix is to derive expressions 
for the Green's function G(R, ui) with R = (m, n), which 
are valid in the uj — > limit. The starting point is Eq. 
@, and the first step is to express 



cos(fc„m) 



~%m — 1 ~~~m 



COS™ k, r 



+^E(-i) j m( r i 2 "v )! ( 2 -M^ 

2 ~[ ]\{m-2])\ 

where [. . .] refers to the integer part of the argument. We 
focus on the half-filled case p = and write Eq. J5J as 
the sum of terms of the form 

g pq = y^ j Coa p (k x )cos g (k y ) UJT " 

k 



where p = m, m — 2, . . . and q = n, n— 2, . . .. We proceed 
by linearizing the dispersion near the node at (7r/2, 7r/2) 
and making the coordinate transformation E 2 = e k + A k , 
tan# = Ak/fik- 



9pq 



271 



d,0 



2p+i ./ 2?r 



sin v 



cos t 
~27" 



smf 



cos t 
~2T 



x [ Ev+a+i ^^ ± E ( cos 6i 3 + sin flfi) 
Jo uj 2 - E 2 

The prefactor is a = N/ (2itvpva) where N — 4 is the 
number of nodes, vp is the Fermi velocity and va is the 
anomalous quasiparticle velocity |VkAk|, and the cutoff 
A is of order Ao. The integrals over E and 8 are easily 
done and 

ffpflM = ^[uF P+q (uj)P° q To+F p+q+1 {u)(p3 q T 3 +Pp q fi)} 

where P^ q are constants given by the angular integra- 
tions, and 



Fa(uj) = / 

Jo 



EdE 



E 2 



The constants P^ q vanish for j — 1,3 when p + q = even 
and vanish for j = when p + q = odd. The first few 
nonzero elements are 



p o 
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Al - ni ~ 2A 2 + 8t 2 
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Only even moments of F a (w) are needed: 



from which we estimate a range of validity 



^1 w 2i A 2(n-i) w 2„ A 2 

F 2n = 22— — + —ln- 



3=0 



2(n - j) 



Since we are interested in the leading order behavior of 
G(R, to) we note that for small w , 



F2n(w) 



1 A 2 
- in - 

2 -w 2 

A 2 " 
~2n~ 



For R = (2m, 2n), the leading order contribution to 
G(R, to) comes from the single term in the expansion 
containing goo- To second order in u>: 



(A/A )2(-" + ") 

\u>\ < Ae 16 (™+") 



For other R, there is no single dominant term in the 
expansion for the Green's function, and the leading order 
behavior comes from the sum over a large number of real 
nondivergent terms. For our purposes, it is sufficient to 
note that when R = (2m + 1, 2n+ 1), the sums take the 
form 



G(R,w) = w<7o(R)fb, 



(16) 



A 2 

G(R, w) = -(-l) n+m ^ In 5 *b + wC (R)fb, (15) and when R = (2m + 1, 2n) or (2m, 2n + 1) 

2 — c<j^ 



where Go(R) is real, and is the sum of several terms. The 
largest term contributing to G(R) is of order 



G(R,w)=G 1 (R)f 1 + G 3 (R)f3, 



(17) 



a\ui\ 



A 



16(m + n) VA 



2(ro+n) 



where Gq(R), Gi(R), and Ga(R) are real constants. 



[1] Ali Yazdani, C. M. Howald, C. P. Lutz, A. Kapitulnik, 

and D. M. Eigler, Phys. Rev. Lett. 83, 176 (1999). 
[2] E.W. Hudson, S.H. Pan, A.K. Gupta, K-W Ng, and J.C. 

Davis, Science 285, 88 (1999) 
[3] S. H. Pan, E. W. Hudson, K. M. Lang, H. Eisaki, S. 

Uchida, J. C. Davis, Nature, 403, 746 (2000). 
[4] T. Cren et ai, Phys. Rev. Lett. 84, 147 (2000). 
[5] S.-H. Pan et al., Nature 413, 282 (2001). 
[6] K.M. Lang, V. Madhavan, J. E. Hoffman, E. W. Hudson, 

H. Eisaki, S. Uchida and J.C. Davis, Nature 415, 412 

(2002). 

[7] C. Howald, P. Fournier, and A. Kapitulnik, Phys. Rev. 

B 64, 1005041 (2001). 
[8] C. Howald, H. Eisaki, N. Kaneko, A. Kapitulnik, 

cond-mat/0201546 
[9] J. E. Hoffman et al, Science 295, 466 (2002). 
[10] J. M. Byers, M. E. Flatte, and D. J. Scalapino, Phys. 

Rev. Lett. 71, 3363 (1993). 
[11] A. V. Balatsky, M. I. Salkola, and A. Rosengren, Phys. 

Rev. B 51 15 547 (1995). 
[12] See M.E. Flatte and J.M. Byers, Solid State Physics 53, 
137-228 (1999). 

[13] Anatoli Polkovnikov, Subir Sachdev and Matthias Vojta, 

Phys. Rev. Lett. 86, 296 (2001). 
[14] J.X. Zhu, C.S. Ting, and C.R. Hu, Phys. Rev. B62, 6027 

(2000). 

[15] I. Martin, A. V. Balatsky, and J. Zaanen, Phys. Rev. Lett 
88, 097003 (2002). 

[16] D. Morr and Stavropoulos, cond-mat/0207243 

[17] Lingyin Zhu, W. A. Atkinson, and P. J. Hirschfeld, 
cond-mat/0208008 

[18] P.J. Hirschfeld and W.A. Atkinson, J. Low Temp. Phys. 



126, 881 (2002). 
[19] N. E. Hussey, Adv. Phys. 51, 1685 (2002). 
[20] A. A. Nersesyan, A. M. Tsvelik, and F. Wenger, Phys. 

Rev. Lett. 72, 2628 (1994). 
[21] T. Senthil and M.P.A. Fisher, Phys. Rev. 560, 6893 

(1999). 

[22] W.A. Atkinson, P.J. Hirschfeld and A.H. MacDonald, 

Phys. Rev. Lett. 85, 3922 (2000). 
[23] K. Ziegler, M.H. Hettler, and P.J. Hirschfeld, Phys. Rev. 

Lett. 77, 3013 (1996). 
[24] C. Pepin and P.A. Lee, Phys. Rev. Lett. 81, 2779 (1998); 

Phys. Rev. B 63, 054502 (2001). 
[25] Claudio Chamon and Christopher Mudry, Phys. Rev. B 

63, 100503-1 (2001). 
[26] M. Fabrizio, L. Dell'Anna, and C. Castellani, Phys. Rev. 

Lett. 88, 076603 (2002); A. Altland, Phys. Rev. B 65, 

104525 (2002). 

[27] W.A. Atkinson, P.J. Hirschfeld, A.H. MacDonald, and 
K. Ziegler, Phys. Rev. Lett. 85, 3926 (2000). 

[28] A.G. Yashenkin, W.A. Atkinson, I.V. Gornyi, P.J. 
Hirschfeld, and D.V. Khveshchenko, Phys. Rev. Lett. 86, 
5982 (2001). 

[29] Inane. Adagideli, Daniel E. Sheehy, and Paul M. Goldbart 

Phys. Rev. B 66 140512 (2002). 
[30] Y. Onishi, Y. Ohashi, Y. Shingaki, and K. Miyake, J. 

Phys. Soc. Jpn. 65, 675 (1996). 
[31] U. Micheluchi, F. Venturini and A. P. Kampf, 

cond-mat/0107621 
[32] Jian-Xin Zhu, D. N. Sheng, and C. S. Ting, Phys. Rev. 

Lett. 85 4944 (2000). 
[331 Y. H. Yang, Y. G . Wang, M. Liu, and D. Y. Xing, 

cond-mat/0211590 (2002). 



12 



[34] R. Joynt, J. Low Temp. Phys. 109, 811 (1997). 

[35] A. V. Balatsky and M. I. Salkola, Phys. Rev. Lett. 76, 
2386 (1996); See also the comment by D. N. Aristov and 
A. G. Yashenkin, Phys. Rev. Lett. 80, 1116 (1998) and 
the reply A. V. Balatsky and M. I. Salkola, Phys. Rev. 
Lett. 80 1117 (1998). 



[36] A. C. Hewson, "The Kondo Problem to Heavy 
Fermions", Cambridge Univ. Press (1993) 

[37] Christopher Mudry, P. W. Brouwer, Akira Furusaki, 
Phys. Rev. B 59 13221 (1999); K. Ziegler, W. A. Atkin- 
son, P. J. Hirschfeld, Phys. Rev. B 64 54512 (2001). 



